Radio frequency spectra of Feshbach molecules in quasi-two dimensional geometries 
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The line shape of radio frequency spectra of tightly bound Feshbach molecules in strong transverse 
confinement can be described by a simple analytic formula that includes final state interactions. By 
direct comparison to experimental data, we clarify the role of effective range corrections to two-body 
bound-state energies in lower dimensions. 
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Recent experimental progress in the creation of quasi- 
two dimensional (2D) atomic Fermi gases Q~HS] has pro- 
vided important insights into the physics of fermionic 
pairing. Pairing phenomena in two dimensions [7] are in- 
triguing because of the presence of a confinement-induced 
bound state for a zero-range interaction potential of ar- 
bitrary strength and sign [5J. This confinement-induced 
bound state significantly affects many-body pairing in 
the BEC-BCS crossover in two dimensions and leads to 
the appearance of a pairing pseudogap phase for weak 
pairing HHn]. 

In the investigation of fermionic pairing, the 
single-particle excitation spectrum [TTHHj probed by 
momentum-resolved [15] or momentum-integrated [21 
HI [B] radio-frequency (rf) spectroscopy plays a crucial 
role. Up to now, theoretical work on rf spectroscopy has 
focused on the strict two-dimensional case whereas ex- 
periments have been conducted in quasi-two dimensions. 
The latter refers to the situation in which tight confine- 
ment along one axis, usually with a harmonic frequency 
cj x , establishes a kinematically two-dimensional gas with 
Ep, ksT <C hjj z . Here, Ep is the Fermi energy of a two- 
dimensional Fermi gas, fcs is Boltzmann's constant and T 
is temperature. However, if the three-dimensional s-wave 
scattering length a s , parameterizing the zero-range inter- 
action potential, is smaller than the axial ground state 
a s < l z = yjh/mu) z , the relative wave function of the 
molecule explores higher transverse harmonic oscillator 
modes. Here, m denotes the atomic mass. 

In this Rapid Communication, we investigate the ef- 
fects of quasi-2D confinement on rf spectra with a par- 
ticular focus on the regime of strong pairing with a s > 
and Eb > fru z , both experimentally and theoretically. 
Specifically, we focus on three different aspects: first, the 
effect of finite range corrections, which are required to 
match experiment and theory in the limit of large binding 
energies, second the bound-bound transitions between 
confinement-induced states, and third the effects of fi- 
nal state interactions on rf spectra, which display signif- 
icant differences between the strict-2D and the quasi-2D 
scenario. Our analytical theoretical results for binding 
energies and line shape functions compare well with the 
experimental data. 

In our experiment, we prepare a 50/50 spin mixture 
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FIG. 1: (Color online) Binding energies of Feshbach molecules 
in a quasi-2D geometry. The bound-state energy including 
effective range corrections according to Eq. (J3j) is shown as 
the solid blue line. The theory for zero-range interaction only 
is the dashed gray line. 

of 40 K atoms in the |F = 9/2, m F = -9/2) = |1) and 
\F = 9/2, rap = -7/2) = |2) states of the hyperfine 
ground state manifold [5] . We achieve quantum degener- 
acy of two-dimensional Fermi gases in an optical lattice 
potential formed by a horizontally propagating, retro- 
reflected laser beam of wavelength A = 1064 nm. The 
trapping frequency along the strongly confined direction 
is lj z = 2tt x 75 kHz, which is calibrated by intensity 
modulation spectroscopy. The radial trapping frequency 
of the two-dimensional gases is oj± = 2tt x 127 Hz, and we 
confine on the order of 2 x 10 3 atoms per two-dimensional 
gas at the center of the trap. Along the axial direction we 
populate approximately 30-40 layers of the optical lattice 
potential with an inhomogeneous peak density distribu- 
tion. Starting from a weakly interacting gas (a s ~ 0) 
at a magnetic field near 209 G, we adiabatically increase 
the interaction strength by lowering the magnetic field 
at a rate of up to 0.25 G/ms to a value near the Fes- 
hbach resonance at 202.1 G. We apply an rf pulse with 
frequency is T f near 47 MHz with a Gaussian amplitude 
envelope with a full width at half maximum of 230 /is 
to transfer a fraction of atoms from the |2) state to the 
|F = 9/2, m F = -5/2) = |3) state. Atoms in the |3) 
state have a two-body s-wave scattering length of 250 



2 



Bohr radii with the |1) state. We turn off the optical lat- 
tice 100 /is after the rf pulse, switch off the magnetic field 
and apply a magnetic field gradient to achieve spatial 
splitting of the three spin components in a Stern-Gerlach 
experiment during time-of-flight. Finally, we detect the 
atom numbers in each of the atomic states by absorption 
imaging. 

We determine the binding energies of the molecules 
by recording the number of atoms transferred into the 
state 1 3) and determine the threshold of the spectrum 
corrected for our spectral resolution of 1.5 kHz [5]. In 
the following, we subtract the (trivial) contribution of the 
atomic Zeeman energy Ez from the rf frequency to ob- 
tain Huj — Ez — hv r f. In Figure[TJwe display the measured 
binding energies and compare with theoretical models. A 
pair of atoms interacting via zero range contact interac- 
tions under quasi-2D confinement has a binding energy 
given by the solution to the equation [3l [8] 
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where g(x) = J °° (l - e— [(1 - e-^)/{2u)} ~ 1/2 

[3]. This prediction is shown as the dashed line in Fig. 
[TJ In the limit of large binding energy Eb > hjj z we 
find significant deviations from the quasi-2D binding 
energy. It has been argued that under tight confine- 
ment and for large binding energies, effective-range 
corrections to the standard zero-range model are of 
relevance [TJ [TjJ [18j l20ti22] . These corrections can be 
viewed as an energy dependence of scattering length 



' eff 



a s (e) 



(2) 



Here e is the energy in the center of mass frame, k — 
y^me/H and r e ff is the effective range. By replacing 
a s — > a s (e) with e = huj z /2 — Eb, the effective range 
corrections to Eq. ([I]) take the form (see also the Sup- 
plemental Material) 
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For open channel dominated Feshbach resonances, the 
effective range is given by [531 HI] 
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1 - 2 



(4) 



where a = 2 7 r(C , 6 m/fi 2 ) 1 / 4 /[r(l/4)] 2 is the mean scat- 
tering length. This implies that, in the vicinity of a 
Feshbach resonance, 40 K atoms have an effective range 
of r c g ~ lOnm (where we used C§ — 3897 in atomic 
units) [5S[ [55] . We plot the binding energy including ef- 
fective range corrections as the solid line Fig. [TJ The 
excellent agreement with our measured binding energies 



demonstrates the importance of effective range correc- 
tions for quasi-2D confinement in the limit of tightly 
bound Feshbach molecules of the |1) and |2) states of 

K atoms. The agreement between and theory and ex- 
periment can presumably be further improved if the fi- 
nite width of the Feshbach resonance is taken into ac- 
count [IH1 [57] . This is expected to give a contribution of 
order —2R* = —2a/s res ~ —a to the effective range [37j . 

Next, we turn our attention to the effects of the quasi- 
2D nature on the rf spectra. First, we study the case 
without final state interactions. This situation is partic- 
ularly simple as it contains only a bound-free transition. 
The rf transition is driven with a Rabi frequency fi/j, de- 
scribed by the operator V = hQ R (\2)(3\ + |3)(2|). The 
transition rate is then given by Fermi's golden rule 

W = Y E I (f\n)\ 2 S(^ + E t - E f ), (5) 

/ 

which involves computation of the Franck-Condon 
overlaps M n {q) = {f\y\i)/{f^R) between the 
normalized bound-state relative wave function 
^ 2 (r) = [nl z g\E B 2 /frw z )]- 1/2 j: n Mz)MO)Ko(KnP), 
with K n = \J m(E 
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Here, </> n (z) is the wave function of the n-th transverse 
harmonic oscillator mode, q is the in-plane relative mo- 
mentum and -E„ q = h 2 q 2 /m + huj z n. After summation 
over final states, we find the resulting line shape to be 
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6(huj - E B 2 - 2huj z j). 



Here, g'{x) denotes the derivative of g(x), 0{x) is the 
Heaviside function. For rf frequencies lo near the binding 
energy, the line shape shows a 8(hu> — E b 2 )/oj 2 scaling 
which is found also in strict 2D [TTJ. In this case, the rf 
spectrum takes the form 



r ( w ) = n R C2D 
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for < huj - E B 2 < 2hw z (8) 



in which we have introduced the two-dimensional contact 
parameter C 2D = 2V2^/l 2 z g'{E B 2 /huj z ) [TTJ US], which 
satisfies the 2D Tan relation. However, as the detuning 
from the binding energy becomes comparable to twice 




FIG. 2: Comparison of the experimental spectra in the molecular regime (a s > 0) with the theoretical line shape functions 
described in the text for three different three-dimensional scattering lengths. We show theory curves for quasi-2D without final 
state interactions (red line), quasi-2D theory with final state interactions (blue line), and the strict 2D expression [see Eq. (13) 
of Ref. [llj ] with 2D binding energies replaced with the quasi two-dimensional binding energies (black line). All curves were 
fitted with the same onset and vertical offset, and they are convoluted with a Gaussian of 2.5 kHz width to account for the 
experimental resolution. Additionally, we have fitted a Gaussian of variable height and width centered at v — to account for 
the atomic peak. The absolute scale of T{v) was least-squares fitted individually for each of the three theory curves. 



the axial oscillation frequency 2uj z , additional local max- 
ima appear in the spectrum, signaling the importance 
of higher axial modes. The tail of T n (uj) at very large 
frequencies Huj 3> EB,fiuj z behaves as 
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using the three-dimensional contact 30J C^d 
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This illustrates the fundamentally three- 
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dimensional character of quasi-2D Feshbach molecules for 
short interatomic distances. Therefore, the strict-2D con- 
tact [11] [29] will only be revealed from rf spectra when 
confinement effects are negligible, i.e. when frequencies 
and binding energies much less than fuv z are considered. 

Finally, we include final state effects. They have to be 
considered if the final state |3) of the rf spin flip is still in- 
teracting with one of the initial states. In our experimen- 
tal situation this will be the interaction between |1) and 
1 3) which is parameterized by the scattering length a^. 
The strict-2D problem of this scenario was studied by 
Langmack et al. |llj and in a three dimensional setting 
by Chin and Julienne [SI] ■ The effect of the final state 
interactions is two-fold: Firstly, a two-body confinement- 
induced bound-state appears and hence we expect the 
spectrum to exhibit always a bound-bound transition in 
addition to the usual continuum [3J. Second, the line 
shape of the continuum spectrum changes as compared 
to the non-interacting case. Applying Eq. [5] (Fermi's 
golden rule), we find that the rf spectrum of a |1)/|2) 
bound-state takes the form 

r(w) = 2nMl 2 R \M\ 2 5(hjjj - E B 2 + E B 3 ) + T c (w) (10) 

Here, M. denotes the bound-bound transition amplitude 



and E B 3 denotes the binding energy in the final state 
channel as a measure of the strength of the final state 
interactions. For the bound-bound transition amplitude, 
we obtain the analytic expression 
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by considering the overlap between the normalized 
bound-state wave functions tp Bl where ij refers to the 
|1)/|2) and |1)/|3) channels (see also the Supplemental 
Material) [31]. Here, we introduced the dimensionless 
binding energies £g = Eg / hw z . Note that this result 
completely includes the effect of the confinement and in 
particular applies to molecules with a characteristic size 
of the order of the transverse confinement or smaller. Re- 
cently, this amplitude was also calculated numerically [B] . 
The strict 2D expression 
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is obtained from Eq. ([IT]) in the limit of small binding 



energies (see also [6]lllj h We find that Eq. (29) provides 
a reasonable estimate for the magnitude of the bound- 
bound transition amplitude even beyond the formal limit 
of its validity E B <C hw z . 

The transition rate into the continuum involves again 
the overlap between the |1)/|2) bound-state and the un- 
bound scattering states of the interacting final state chan- 
nel. With the analytic expression of Ref. [S] for these 
eigenstates ^|(r) in the presence of short range inter- 
actions between |1)/|3) atoms, we obtain the expression 
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(see also the Supplemental Material) 

M n (q) = f d\ [< 3 g (r)]*V#(r) 



/ 12 (£„ q ) 
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which contains the scattering amplitude between atoms 
in states \i) and \j) at relative energy E [5]: 



,)-g{-E/Hu z )' 



(14) 



After summation over final states in Fermi's golden rule, 
the transition rate into the continuum takes the form 



r c ( w ) = roM-FH, 



(15) 



where effect of final state interactions on the continuum 
spectrum is described by the factor 



Eh 2 ) 



(16) 



When considering dissociation into the lowest axial 
mode at sufficiently low energies E <C Sw z , the scattering 
amplitudes f* 3 (E) may be approximated as 



f\E) 



in 



2irl z /oij + ln(-Bhuj z /{irE)) ' 



(17) 



where B ~ 0.905 J3J IS] - In this approximation, J-(oj) 
becomes 



ln[e 12 /(fiw - £#)] 



ln[c 13 /(fiw - E%)} + : 



(18) 



with e»j = (B/Tr)huj z e^ rl ^ aij . The energy scales ej 



appearing in Eq. ( 18 1 are analogous to the scattering 
length a.ij in three dimensions in the sense that they are 
the single parameter describing low energy collisions be- 
tween atoms in the hyperfine states \i) and \j). How- 
ever, note that e.y is different from the binding energy 
for deeply bound molecules under quasi-2D confinc- 
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ment. We emphasize that Eq. ( 18 1 holds for any interac- 



tion strength at low energies (i.e. as long as E <C Hu) z ). 
This line shape is in general different from what would 
be obtained by using the strict two-dimensional formula 
of Refs. |H |TT] and replacing the 2D binding energies 
with the actual quasi two-dimensional binding energies 

£,12,13 
B ■ 

In Figure [2] we display our measured rf spectra and 
compare with three different theoretical models: strict- 
2D and quasi-2D with no final state interactions, and 
quasi-2D with final state interactions. Generally, we find 
that the quasi-2D theories are more strongly peaked near 



the binding energy which appears to fit the data bet- 
ter. This difference in line-shape is particularly strik- 
ing when the initial state is a tightly bound molecule or 
if the interactions in the final state are strongly repul- 
sive (i.e. ai3 3> l z ). In contrast, for our experimental 
data for 40 K with weak repulsive final state interactions 
(W a i3 ~ 4.61 ± 0.05 over the range magnetic fields con- 
sidered here), a quasi-2D line shape that does not take 
into account final state interactions still provides a good 
fit, whereas the strict-2D theory suggests a lesser-peaked 
line shape. The recent Duke experiment [5], on the other 
hand, dissociates a weakly bound molecule into an out- 
going channel with weak attractive interactions. In their 
case, E^ (E^) approximately coincide with ei 2 (£13), ex- 
plaining why their strict-2D approach gives a good result 
for the line shape of a quasi-2D system. 

In order to gain deeper insight into the more compli- 
cated and fundamentally important many-particle prop- 
erties of confined quantum gases |3"2TI3"4"] , a firm un- 
derstanding of the interplay of confinement and few par- 
ticle physics is crucial. In this Rapid Communication we 
provided analytic results for rf spectra of dilute quasi-2D 
paired gases in the two-body limit and compared those 
results to experimental data. Additionally, we demon- 
strated how effective range corrections contribute to the 
binding energy of confined quantum gases by explicitly 
comparing theory to experimental data. Beyond quasi- 
2D quantum gases, we believe our results can be gen- 
eralized to quasi-one dimensional paired Fermi gases as 
experimentally studied in Refs. [35 , 36 . 
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Supplemental Material for Radio frequency spectra 
of Feshbach molecules in quasi-two dimensional 
geometries 

Binding energy and effective range corrections 
under confinement 



Theoretically, two interacting atoms in hyperfine states 
|1) and 1 2) under quasi-2D harmonic confinement are de- 
scribed by the Hamiltonian 

h 2 dl , i 

+ V H (zi) 



H = 



E 



2m 



+ V- nt (ri-r 2 ), (19) 



where the interaction Vi n t(r) is the interaction potential 
and Vh(z) = \mio\z 2 . After separating into the relative 
coordinate r = ri — r 2 and the center of mass coordinate 
R = (ri + r 2 )/2 one has 
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Muj 2 Z 2 



(20) 
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(22) 



with reduced mass \x — m/2 and total mass M = 2m. 
For a short ranged potential Vi n t(j) with range tq -C l z , 
we can find an approximate solution to the two-body 
problem by first noting that the ansatz pQ 

fe(r)«G £ (r,0) (23) 

where Ge(?i, r 2 ) is the Green's function of the Hamilto- 
nian H re i. G E (r, 0) satisfies the Schrodinger equation for 
the relative coordinate Hamiltonian away from the origin 
r = 0. Close to the origin (tq < r -C l z ) we match this 
outer solution to logarithmic derivative of the asymptotic 
free space scattering solution of the interaction potential 
V mt {v) 

d r (rtp E (r)) 



rip E (r) 



(24) 



r=0 



y/2/j,E/H and S E is the 



k E cot(5g 

where E = E + huj z /2, ki 
s-wave phase shift at energy E. With the usual effective 
range expansion k cot [6) = —l/a s +k 2 r c ^/2- 
expansion of G E (r, 0) from Refs. 013] 

g(~E/hu z ) 



and the 



G E (r,Q) 



L 



(25) 



we arrive at eq. 3 of the main text. 



Effect of confinement and final state interactions on 
radio frequency spectra in quasi-2D 

Simple considerations: 2D Limit 

It is instructive to study the problem of molecule disso- 
ciation in a strict two-dimensional setting jH [S] . Such a 



description should be valid whenever the binding energy 
E B 2 (E B 3 ) of the molecule in initial (final) state chan- 
nels and the energy of the RF-photon is much less than 
the spacing of transverse harmonic oscillator modes uj z . 
In general the RF spectrum T(oj) of molecules can be 
calculated from Fermi's Golden rule [B] 

T 2D (u) = 2nJ2\(^f)\ 2 6(u J - E B 2 E f ) (26) 
/ 



The summation over final states includes a bound-state 
at energy Ef — —E B 3 and a continuum of unbound atoms 
with Ef = k 2 and relative momentum k. For brevity, we 
use units where h — m = ttji = 1. The matrix element 
for the bound-bound transition amplitude M2D is given 
by the overlap integral 



M 



2D 



d 2 r 



[^ 2 (r) 



r) 



(27) 



where ip B 2 (r) [V'b^ 1 ")] is the molecule wavefunction in the 
initial [final] state, respectively. At bound-state energy 
E B = K 2 j , the normalized wavefunction for short-range 
attractive s-wave interactions in 2D is 



V^(r) 



KijKoiKijp)/^ 



which gives 



M 



2D 



^EjE£\og(E B 2 /E B *) 



~ E B 3 
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(29) 



We note that this result makes sense from a qualitative 
point of view: When interactions in final and initial state 
are equal, one has M2D = 1 (as all weight is concentrated 
in the bound-bound transition), whereas for vanishing 
final state interactions (i.e. E B 3 = or E B 3 = 00), M2D 
tends to zero. 

The bound-free continuum can be calculated from the 
matrix element between a s-wave scattering state |^ + ^) 
and IV'b)- Such a state satisfies a Lippmann-Schwinger 
equation 



i^ +) > = \4 0) ) + 



1 

k 2 — Hq + irj 



(+)\ 



(30) 



and is guaranteed to have the same normalization as 
the solutions to the free Hamiltonian 



For a 2D 



short range potential the Lippmann-Schwinger equation 
is solved in the s-wave channel by [H [7] 

d + )r^-M¥l_^i_ fl3{k 2 )H ii) {kp) (31) 



2tt 2tt4" 



with the 2D scattering amplitude for the final state chan- 
nel 



r{E) 



4tt 



in + \og(E B 3 /E) ■ 
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Here we picked the normalization of the scattering states 
such that 



<.Mu\ * ./,(+) M - J_ S(k - k') (33) 

2ltk 



4 +) (r) ^(r) 



In this normalization, the bound-free transition ampli- 
tude M.2D (k) becomes 



M 2D (k) = / d\ [^ 2 (r)]*4 +) (r) 



k 2 + E 12 



1 - 



f 12 (E) 



(34) 



Evaluating the sum over final momentum states we ob- 
tain a result for the bound-free spectrum part of the spec- 
trum T 2D {uj): 

Tl D {uj) = 2nJ d 2 k \M 2D {k)\ 2 6(Lj - Eg - k 2 ) 

= r^H^M (35) 

where T 2D (lu) = 6(uj - E 12 ) 2 ^^- is the RF-spectrum 
without final state interactions and 



f 13 (u-Eh 2 ) 



(36) 



is a correction term accounting for final state interac- 
tions. Note that the complete spectrum 

r 2 » = \M 2D \ 2 5(E™ El 2 + Tl D {u)T 2D {u) 

satisfies the sum-rule 



doj T 2D (lj) = 2ir. 



(37) 



Effect of closed channels 

We will now turn our attention to the problem of a 
three-dimensional gas confined into two dimensions by 
tight one-dimensional harmonic confinement. In such 
a potential the relative and center of mass coordinates 
separate, so we only need to consider the Hamiltonian 
for the relative coordinate. In order to calculate the 
bound-bound and bound-free transition rates, we first 
need to know the normalized wavefunction of the two- 
body bound state in this geometry. The Green's function 
at energy E in the harmonic oscillator potential is given 

by m 

where q 2 = E — v and the transverse harmonic oscillator 
states have 0„(O) = (y- l)ll/[V^.(2irl 2 z ) 1 / i }. We can use 



the explicit expression for the bound-state wave function 
ijjg (r) = NijG_ E ij (r) in terms of G.e(r, 0) to obtain the 
normalization factor 



-in 



(38) 



where Eg is the two-body binding energy, satisfying 
g(Eg /huj z ) = l z /a,ij. The normalized wavefunction of 
the two-body bound-state is thus given by 



^nl z g'(E^/uj z ) v 



With this wavefunction it is straightforward to calculate 
the bound-bound transition amplitude M. 



M 



l^(0)| : 



,\og{\q u \ 2 /\q' v \ 2 ) 



g(& 2 ) -g(ffi) 
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K\ 2 



(39) 



where £^ = Eg / ftw z . We note that for equal interactions 
in initial and final states (i.e. ^3 = £,]^), M. = 1. Also 
note that for small £]j 2 , £]j 3 <C 1, eq. (39) reduces to the 



2D result eq. (29). Analogous to the simpler 2D limit 



we can obtain the contribution of the bound-free contin- 
uum from the overlap between scattering state and bound 
state wavefunction. For the final scattering states we use 
the basis of retarded scattering states given described in 
Ref. 



1 

2^ 



(j) u {z)J a {qp) 



^G E (r,0) (40) 



For the precise expressions of the coefficients A' v see 
Ref. [2J. Here the prime indicates that we consider the 
final state channel. By parity symmetry, the overlap ma- 
trix element AA v (q) vanishes for odd v. For even v one 
obtains 

MM) = I d\4> 12 {r)^+ q \r) 



MO) 



f 13 (E) 



f 12 (E) 



(E-E$Wnl,g'(E)*/u> a ) 

E=q 2j r v 

The transition rate into the continuum is given by 
r c (w) = 2TrJ2\Mv(q)\ 2 8(u-E 1 I ?-q 2 - V Lj z ) 



= 2t: 2 Y,\M(Ju J -^ z -E} 2 )\ 2 



'2n 



E 



(2J-1)! 



■6{uj - E 12 - 2j) 



/ 13 (u, 



Eh 2 ) 



/12( W 



Eh 2 ) 



Note that this result is again of the form T c (uj) — 
ro(w)J r (u;), where ro(w) is spectrum without final state 
interactions and the factor 



after exchanging integration and differentiation. We can 
then numerically compute 



Eh 2 ) 



Eh 2 ) 



(41) 



We finally note that it is convenient to use a closed form 
expression for the sum 



E 



(2j-l)!! , 
(2j)!! 



(x - 2j) 



2 r([x/2\ +3/2) 



(42) 



as given in Ref . ;2] . With this expression one can directly 
confirm that the tails of RF-spectrum without final state 
interactions obeys the 3D Tan relations. 



W(e 1 ,e 2 )= / g'(z)dz 



(45) 



where C is a contour in the lower complex halfplane (e.g. 
we used straight lines along e\ — » —i — » —62). 



Numerical evaluation of g(e) 

Since the expression for <?(e) via the summation of Ref. 
[2] converges very slowly we will give a numerically more 
efficient scheme to evaluate differences of the form 

W(e 1 ,e 2 )=g(e 1 )-g(-e 2 ) (43) 

where £1,62 > 0- Using the integral representation for 
g'(e) of Ref. [3] one obtains 
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